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ABSTRACT 


A computer program based on the theoretical work of Gyftopoulos, 
Steiner, and Levine on bimetallic systems and using a modified version of 


This program, WFGSL, accepts the operating conditions and the physical 
parameters pertinent to the substrate and adsorbate and outputs the field - 
free work function, electron current (Richardson equation), ion current 
(Saha equation), and fractional substrate coverage by the adsorbate. A brief 
description of the theory is presented together with a program description 
and listing. An application of the program to a bimetallic system of cesium 


(adsorbate) and rhenium (substrate) is also described 
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I. INTRODUCTION 


In thermionic energy conversion, the work function of an emitter in 
contact with cesium vapor is a controlling parameter determining the rate of 
electron emission and thus converter efficiency. An understanding of the 
adsorbate- substrate interaction is essential to the development of more effi- 
cient thermionic converters. The theoretical work of Gyftopoulos, Steiner, 
and Levine (Refs. 1, 2) on the general bimetallic system evaluates the work 
function of this system as a function of its physical parameters. 

This theory has been incorporated in a program using a modified ver- 
sion of Wilkins' SIMCON subroutine SURFAS (Ref. 3). The program serves 
as a tool in evaluating possible changes in a converter due to changes in the 
work function of the emitting surface with time. Although cesium is the only 
adsorbate considered here, the computer program is general in nature and 
can be used for other bimetallic systems. The program accepts the operat- 
ing conditions and the physical parameters pertinent to the substrate and 
adsorbate and outputs the field-free work function, electron current, ion 
current, and fractional substrate coverage by the adsorbate. The code 
differs from Wilkins' program in format usage, input parameters, and out- 
put information; the printed plus graphical output facilitates the comparison 
of theoretical and experimental work. 


II. THEORY 

The work function of a pure metal can be identified with the neutral 
electronegativity of valence orbitals in the surface atoms. Using this ap- 
proach, the work function for a bimetallic surface can be viewed as a per- 
turbation from the bare work function. The perturbation arises from 
substrate-adsorbate (s-f) and adsorbate-adsorbate (f-f) interactions (Fig. 1). 
In an s-f bond, there is a charge transfer due to a difference in the 
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electronegativities of the atoms. This charge transfer results in a dipole 
double layer whose electrostatic potential is superimposed on the electro- 
negativity of the substrate atoms. The perturbation which results is taken 
to be proportional to the electronic charge transferred, F. The f-f bond 
resulting from charge overlap also produces an electrostatic potential, again 
affecting the bare substrate electronegativity. This perturbation is assumed 
proportional to the amount of charge overlap, Q. Thus, our basic relation 
(Refs. 1, 2) is 


<b = <i> + cQ + bF, 

s 


( 1 ) 


where 4> is the perturbed neutral electronegativity, cji s I s the bare substrate 
work function (unperturbed neutral electronegativity), and c and b are con- 
stants of proportionality. Observe that 4> = cja s at 0 = 0 (0 is the fraction of 
substrate coverage by the adsorbate) since Q and F are zero by definition. 

It is found experimentally that with a one-monolayer coverage (0 = 1), 
the work function is approximately equal to the adsorbate work function, cj^.: 

<(j = <t>£, Q = 1 • 


Equation (2) is therefore a boundary condition for (1). 

The quantities cQ and bF may be derived, and result in the following 
equations: 

cQ = -(<j> g - 4> f )M, (3) 


where M is a Morse function defined by 


M = 2 exp 


+ad Q (l 


,- 1 / 2 , 


]- 


exp 


^+2ad Q (l - 


,- 1 / 2 , 


(4) 


where ad^ = 2. 97. 
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For bF it is found (Ref. 1) that 


b 


. - 6 


-0.905 (10 ) V f (R cos (3) cr f 0 


(' • ?)( 


1 + 9ntr^ 2 0 3//2 


) 


e V/coulumb, 


(5) 


(<|> - <j> f )(l - M) 



( 6 ) 


where is the number of bonding orbitals per adatom, R and p are defined 
in Fig. 2, cr is the adsorbate density, D Q is the energy of the purely co- 
valent s-f bond, and is an adjustable parameter, which may be approxi- 
mated by (Ref. 1) Dj = 1.3 (4> g + <t> f ) (-28.75/R), where R is in A. Note 
that (2) is satisfied since F(l) = 0. 


III. PROGRAM 

The program logic determining the cesiated work function cj) and the 
fractional coverage 0 as functions of Tg and consists of two nested 

iteration loops. The inner loop uses Newton's method to solve for F in (6), 
assuming an initial value of 0. 5 for 0 and 0. 5 for F. After solving for F, 
the arrival rate is set equal to the desorption rate of adsorbate atoms. 

Since the latter quantity is a function of both cj> and 0, it will not equal the 
arrival rate for the initial guess of 0. Thus 0 is changed accordingly to 
bring the arrival and desorption rates into equality. Having adjusted 0 
according to the difference between the arrival and desorption rates, it is 
now necessary to recalculate F to determine the new value of tj>. Figure 3 
presents a flowchart of the SURFAS subroutine. Table 1 presents the input 
parameters. Besides the output listed in the introduction, WFGSL uses the 
CalComp plotter to generate a graph of cj> versus T g/ T R - A listing of 
WFGSL is given in Appendix A, along with a sample data deck and the corre- 
sponding output (Appendix B). 
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IV. RESULTS 


The results appearing in Figs. 4-8 are for a rhenium substrate and 
cesium adsorbate. The physical constants that were used appear in the data 
deck of Appendix A. The curves in Figs. 4-7 show the dependence of the 
work function cj> on the parameters cM » c , & , D-, ( bar e work function, 
substrate density, adsorbate density, and adjustable parameter, respectively). 

In Fig. 4, <j> is plotted as a function of the ratio between substrate and 
adsorbate temperature for various values of the bare work functions. As is 
expected, a change in the bare work function results in a vertical displace- 
ment in the curves. Figure 5 shows the effect of varying the substrate 
density cr g while keeping the adsorbate density fixed. The effect is seen to 
be most pronounced for smaller values of T g/ T R > where the fractional cover- 
age and thus the substrate -ad sor bate interactions are greatest. In Fig. 6, 
the substrate density is held constant while the adsorbate density is varied. 

The curves show the strong dependence of cj> on the adsorbate; with increas- 
ing adsorbate density the work function is lowered significantly even at high 
substrate-adsorbate temperature ratios. Figure 7 shows that cesiated work 
functions shift by a constant amount over the entire range of temperature 
ratios as the adjustable parameter is varied. In Fig, 8 an experimental 
curve (Ref. 4) is shown together with a theoretical curve generated by 
WFGSL. The agreement is seen to be excellent. 


V. CONCLUSION 

A computer program, WFGSL, based on Gyftopoulos, Steiner, and 
Levine's work on bimetallic systems and using a modified version of Wilkins' 
SIMCON subroutine SURFAS has been developed for use on the UNIVAC 1108. 
The program determines the field-free cesiated work function, electron cur- 
rent, ion current, and fractional substrate coverage by the adsorbate. The 
indicated relation between cf>, the calculated work function, and the various 
physical input parameters is found to be consistent with expectations; 
agreement between experimental work functions and theoretical work func- 
tions is shown to be good. 
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Table 1. Input parameters 


Card 

Columns 

Parameter 

Units 

Format 

1 

1-5 

Number of independent cases > 1 

— 

15 


6-80 

Blank 



Parameter, substrate properties 

2 

1-8 

Atomic weight 

amu 

F8. 3 


9-16 

Bare work function 

e V 

i— ' O o 

r o . d 


17-24 

Full coverage work function 

eV 

F8. 3 


25-32 

Heat of sublimation 

eV 

F8. 3 


33-40 

Metallic radius 

0 

A 

F8. 3 


41-48 

Bond strength 

— 

F8. 3 


49-56 

Substrate atom density 

A" 2 

F8. 3 


57-64 

Adsorption site density 

A" 2 

F8. 3 


65-72 

D i 

eV 

F8. 3 


73-80 

Blank 


F8. 3 

Parameter, vapor properties 

3 

1-8 

Atomic weight 

amu 

F8. 3 


9-16 

Heat of sublimation 

eV 

F8. 3 


17-24 

Metallic radius 

0 

A 

F8. 3 


25-32 

Bond strength 

— 

F8. 3 


33-40 

Ionization potential 

eV 

F8. 3 


41-48 

Electron-neutral cross section 

A 2 

F8. 3 


49-56 

Ion-neutral cross section 

A 2 

F8. 3 


57-80 

Blank 


F8. 3 
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Table 1 (contd) 


Card 

Columns 

Parameter 

Units 

Format 

4 

1-5 

Number of substrate temperatures to be read 

— 

15 


6-80 

Blank 



5 

1-10 

Substrate temperature T g 

K 

F10. 5 


11-15 

Number of adsorbate temperatures to be read 
in for the given substrate temperature 
(NURT) 


15 


16-80 

Blank 




Parameter, NURT 

6 

1-10 

Adsorbate temperature T^ 

K 

F10. 5 


11-80 

Blank 





Fig. 1. Bonding interactions 




Fig. 2. Maximum ligancy on the 1 1 1 -plane 
of tungsten 
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Fig. 8. Rhenium work function 
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APPENDIX A 


PROGRAM LISTING 
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DIMENSION RATI 10 ) ,WF ( 10 ) 

' ' COMMON AC- , HGVP.GVS'G , V I G'» VALGVCSE'GYCSI G " 

REAL JO, 10 
V ALG- 1 • 0 
CALL PLOTS 

CALL PLOT(0.0,C.C,-3) 

READ ( 5 » ICO > NUCS 

■”'5oT0 ■ NUCSNVr.NUCS' 

R e AD t 5 » 10 1 ) AS:FOS;FCS, HS , RS , SS » SDS s AD 0 - » C$ 

" "r"f‘a5~ c by re it " a'g ~,h gt r g" ; s~gvv tg , cs eg vc stg 

RE AD I 5 * 100 ) NUSS 

— wTTTTFTSTITcn 

Wft I T E ( 6 » 1 0 3 ) AS , FOS » FC s , HS » RS , SS , SDS , ADS , CS 

'"'wnr rsvmy 

WRITE ( 6, 103 ) AG ,HG»RG*SG, VIG,CSEG,CS IG 

'■"w^rreiyvrr?) 

DO 10 NUSBCN= 1 , NUSS 
R'E A P ( D , 1 0 2 ) T S , NUR T 
I=NURT+1 

- j = r+ i 

\ W F ( I ) = 2 • 0 

WF"( J )"='C Y2S 

FAT( I ) =2.0 
RATI J)=0.25 

I F ( NUSBCN . NE._1 ) GO TO 3 

'""CA’lU'PLOT ('12".0»"6.6"*-3') 

CALL AX IS(0.0,C.0*5HTS/TR, -3,10.0, 0.0, RATI I ) ,RAT(J) ) 

'"CALL" AX IS (0.rf,'0. O', 1 ?'HWGr 7 ~F 0 NC Tf ON » T 3 V 1 0~ • oY 9~C Yo‘, wPTl" T,'w'F' ( ji") 

3 CONTINUE 

PO"! NURTCN=1,N'JRT 
READ! 5 ,102) TP 

" "p G'= 2".' 4~E 1 V 3 * E XP ( 8 9 1 b’ .' 67 TP ] 7 S OR T 7 T Y ) 

CALL SURF AS ( AS , F OS , FCS ,r S , RS , SS , SDS , ADS , C S , T S , PG , Wr SO , JO , 1 0 * DS , 

' XNOCONV) 

RAT I 0 = T S/TR 

IF ( NCCG.N V . NE . 1 ) GO TO A 
WRITE (6,107) TS»TR 

""'WF' (iNUR'T CN ) ''=27 o'” 

_ _ p AT JA IJ p L- N > = R a t i_ c 

' N0C0NV""=" C 

GO TO 9 

4 CONTINUE 

RAT ( NURTCN ) =RAT 10 

wf i n' Or' f c'n ') v w f' sc" 

I F (NURTCN. NE. 1 ) GO TO 3 


1 z 
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s 

-ta-RIJ-E 1-6 j-10-S ) mS. > J ->J3 AJ. 1X1 > lOfc SO » -LU j.. I LU US 

GO TO 9 
CONTINUE 





9 

WR I T E ( 6 * 1 06 ) TR , RATIO, WFSO . JO. 1 0 . DS 
CONTINUE 





] c 

CALL LINE! RAT >WF ,NURT» 1 .- 1 .NUSbCN ) 
CONTINUE 





ICC 

CALL PLOT) 10.0.0.0.999) 
FORMAT (15) 





101 

102 

FORMAT! 10F8.3) 
FORMAT ( FIG. 5 . I 5 ) 





103 

104 

FORMAT! 1H0.9F12.3) 

FORMAT! lHC.lbH GAP SPACING =»F5.1) 





120 

FORMAT (108H0 AT WT B W FUN C 

lAT'IUS BOND STR SUB DEN AOS DEN 

V; FUN SUR H 

01 ENERGY) 

EAT 

R 


121 

FORMAT I 5 4 H 0 AT WT SUB HEAT 

IN POT E-N XSEC I — M XSEC) 

RADIUS BONO 

STR 

10 


122 

FORMAT ( 8 4 i i 0 SUB TEMP RES TEMP 

1 JO 10 FRAC COV) 

TS/TR WORK 

F UN 



105 

106 

FORMAT ( 1HG . 5F 12 . 3 , F 1 2. 6 , F 12. 4 ) 
FORMA T ( 1 3 X . 4 F 1 2 . 3 » F 1 2 • 6 , F 1 2 . 4 ) 





1 0 7 

FORMAT ( 1CX *3CHSURFa'S fS UnCONV £ RGEO FOR 
3 5.2! 

T S = * F8 • * 1 OH 

AND TR r 

9 F 


-FOR , 

END 

IN SURF AS 






SUSP CUT" I NE SURF AS ( «S ."EOS .FOG.HS.KS. S"s , F.S'S 
1S.N0CC.NV> 

.ASF .C.TS.PG. FS 

.EES.EIS 

»C 


C 5. UR FAS SURFACE PHYSICS SUBROUTINE 

COYMCN AG * HG » RG » SG > V 1 6 * V AlG> CSE o >CSIG 

cc 

c 

EVALUATE ADSORPTION SYSTEM PARAMETERS 





cc 

S'a'G = 2.0 





SDS = 1C0.C*SSS 
ADS = 1 00 . C *ASS 

R = RG + RS 

COSR = SORT ( 1 .C - 5C.G/(SDS*R**2.0) ) 

ALF = KG*~*3 + i.65*RS**2 

Cl = 0.‘>05«R*ADS*CCS3/( l.C+ALF/R**3> 

C? = 9.0E-3*ALF+ADS** 1 . 5 

<- -> - r 

S = 2 . r / ! SG/S5 + S S/S6) 
DO = $QRT(HS*HG)*S 

L L 

A S - 0 • -> 6 £ + A * p G / S ( ». R T ( A G * T S ) 





C C 

BEGIN SURFACE COVtKAGt ITERATION 
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N-C-5-=--0 

CS = 0.5 

D£L_CS__=__Q._5_ 

10 CONTINUE 

TlC.S-.5_ HC.S._±^I 

I F ( NCS.GT . 35 ) GO TO 45 

G__=_ _1 ._0 _ _ _3_ . Q* CS * C S_ _+_ 2.0 *C S*C5* CS 

OG = -6 .0*CS* ( 1.0-CS) 


C BEGIN FRACTIONAL CHARGE TRANSFER ITERATION 

CC 

NF = 0 

F r_Q ;A 

12 CONTINUE 

NF_ = NF +1__ 

IF ( NF . GfVl 5 ) GO “TO' 45'" 

FOF = DO*F/SQRT( l.Ci-F*»2) + C3*F - (F3S-FOG)*G 

DFOF = DO/( 1.0-F**2>**1.5 + C3 
FN = F - FOF / DFOF _ 

IF" ( F'nVl E. 6V6T F N' = ' 0.5*F 

IF(FN.GE.l.O) FN = C.5M1.0+F) 
f F"("a"BS"( FN-F ) '. LEVoVoiT Go"? 0 " 15''" ' 

F = FN 

GO TO 12 


C END FRACTIONAL CHARGE TRANSFER ITERATION 

CC 

15 fa's = DO*SGRT ( 1 ,0-F**2 )' - 0.5*C 3*F** 2 + F* ! FOS-FOG ) *G 

OF AS = F*( FOS-FOG ) »DG 

FS = FOG + ( FOS-FOG ) *G - F*C1 *CS / ( 1 . 0 + C 2*CS** 1 . 5 ) 

FREQ 1 = 2.50E-6*SGRT(FAS_1*SGRT(J AG+4 . 0*A_S )_/ (_2._0*AG*AS ) )/<R*COSo> 

U = " O’." 5* CS* OF AS /FAS 
FR EG_ = FREQ1*EXP ( U ) 

ENT ="'e"xp'(' 6".'^*('CS 6'-CS'r+SQR'f'('crs')'/ ( "l'.'o'-s'oR T ( CS)')') ) / ( SGRT'i l76-CS7 

1*( l.O-SGRT(CS) ) ) 

ERS = 1.0E+14*FRE&*SWG*ADS*CS*ENT*EXP(-11600.0*FAS/TS) 

I F | ABSJ ( ERS-ARS ) /ARS ) -0 .00 5 ) 40 » 20 »20 

20 IF ( E' R S - A R S)'25','40 ,'E 0 

25 DELCS =_ DELCS/2.C 

CS’=‘ CS + DELCS ' ' ‘ ' 

GO TO 10 

SO DELCS = DELCS/2.C 
CS = CS - DELCS 

go 'to' ’ id' ’ 

CC 

C E N r D "S UR FA C~E~ ' CO VE R A G tT ' "I f E R AT I ON " ‘ " 

CC 
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C _ COMPUTE .ELEC.T RON ANO__IO,N. _EM I.SSION.CU RRc.NL.P.FN.SI.T IES 

CC 

40 EES 12_0.0*TS_**2*EXP( - 1 1 6CC_.G*FS/TS_) 

p"f s~ ~=~ ~2 ~.S 9~8 E + 1 * P G~* Efx P(1 160C'.d*-rFS-VIG) /TS ) /SORT ( AG* T S/ 1 1600 . 0 ) 

GO TO 50 

45 CONTINUE 
NOCONV = 1 
50" CC NT I NU E~ 

RETURN __ 

-MAP 

LIB LIB*PLOTS 
— XOT 

- 

186.2 4.8 1.95 8.07 1.35 2.62 0.101 0.029 0.72 

"T 3~3 • 0 0.7 3 2 lb 6 T.O 3.89 40 0. 5'70.0 

5 ; „ 

2043. 3 

486.4 
475 "Y 

464.3 

1940." 4 

485.0 

473.2 
461.9 

45Y.2 

184C.0 5 

y-y-y 

497.3 

484.2 

471.. 8 

4 6 O’ 

1740.0 5 
5‘fiV8 

497.1 

433.3 

470.3 
45 7". 9" 

1660.0 7 

535"5— 

518.8 

503.0 

488.2 

4 7 4” 3- 

461.1 

-44-8.6 

-FIN 


N I F — 
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APPENDIX B 


PROGRAM OUTPUT 
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iXif 







AT WT 

B l* FUN 

C ■*' FUN 

SUB HEAT 

RADIUS 

BOND ST 6 

SUP, OF N ADS OF N 01 F N F ROY 

i 36. 20 3 

4. $00 

1.950 

2 .0 70 

1.3 50 

2.520 

.101 ,029 .720 

AT ut 

SUB HEAT 

Radius 

ECN 3 SIR 

ION POT 

E-N XS ITC 

J-N XS F C 

1 33. oao 

. 730 

2.36 0 

1 .3 00 

3.8 90 

400.000 

570.000 

SUE TEMP 

RES TEMP 

TS/TR 

WCFk F’Jn 

JO 

10 

F P A f CD V 

2343. GQG 

436. 430 

4.20C 

4.139 

.031 

2.714999 

.1118 


4 75. ICO 

4 . 3CG 

4 .260 

.016 

3.534752 

.0520 


454.300 

.4.400 

4.378 

.005. 

4 * 5 p 3 p.l 9 

. 0 7 .3 7 

1940. COO 

465.000 

4 • C 0 D 

2 • 0 3 5 

.037 

.625182 

.15 2 3 


473. 2 jO 

4.10 c 

4.317 

.017 

.832868 

.1313 


4E1.2C0 

4 . 2 DO 

4.146 

.00 8 

1.231156 

.110 4 


451.230 

4. 300 

4.279 

.003 

1 .728.322 

.0893 

1B40.D0C 

5 11.1 00 

3 .600 

3.360 

.257 

.C58C77 

.235 4 


497, 330 

3.70 0 

3.432 

.119 

.078551 

.2153 


464.200 

3.6 00 

3.611 

f; 5 T 

.110640 

.1 94 .3 


471.330 

3. 30 C 

3.752 

.022 

.157982 

.1733 


4 SC. 000 

4 .0 00 

3 • h 8 8 

. CC9 

.247249 

.1519 

1740. coo 

TiiTsoo" 

3 • ^ u C 

3.103 

.365 

.009419 

.2773 


497. ICG 

3 . 50 c 

3.22 7 

.165 

.012614 

.257 3 


433. 330 

3 . 6 QC 

3.351 

.072 

.017518 

. 2358 


4 7C. 300 

3 . 7 CD 

3 • 4 7E 

. n 3 1 

. P 24 E 1 1 

. ? 1 E 3 


457.930 

3.30 0 

3 .£08 

.013 

.035023 

.19 4 3 

lesc.oco 

535. 50G 

3.10C 

2.756 

1.430 

.001355 

.3386 


513.303 

3.20C 

2.358 

. 702 

.001642 

. .313,1 


5 C 3 • Li C u 

3.3GC 

2.581 

.298 

.002294 

.299 3 


438.230 

3.400 

3.238 

. 1 31 

• Q G 3 G 3 0 

. .2.791 .. 


4 74.30C 

3 . 5 CD 

3.218 

.05 7 

.004245 

.2588 


451. 130 

3 . 0 3 C 

3.343 

__ ..3 24. 

. . . .005035 

.2330 


448. COO 

3 . 7CC 

3.4 73 

.CIO 

.008646 

.2168 





